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Abstract 

We study the dynamical magnetization process in the ordered ground-state phase of the trans- 
verse Ising model under sweeps of magnetic field with constant velocities. In the case of very slow 
sweeps and for small systems studied previously (Phys. Rev. B 56, 11761 (1997)), non-adiabatic 
transitions at avoided level-crossing points give the dominant contribution to the shape of mag- 
netization process. In contrast, in the ordered phase of this model and for fast sweeps, we find 
significant, size-independent jumps in the magnetization process. We study this phenomenon in 
analogy to the spinodal decomposition in classical ordered state and investigate its properties and 
its dependence on the system parameters. An attempt to understand the magnetization dynamics 
under field sweep in terms of the energy-level structure is made. We discuss a microscopic mecha- 
nism of magnetization dynamics from a viewpoint of local cluster flips and show that this provides a 
picture that explains the size independence. The magnetization dynamics in the fast-sweep regime 
is studied by perturbation theory and we introduce a perturbation scheme based on interacting 
Landau-Zener type processes to describe the local cluster flip dynamics. 

PACS numbers: 75.10.Jm, 75.50.Xx, 75.45. lj 
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I. INTRODUCTION 



It is well-known that non-adiabatic transitions among adiabatic eigenstates take place 
when an external field is swept with finite velocity^. In particular, at avoided level-crossing 
points strong non-adiabatic transitions occur, causing a step-wise magnetization process^. 

In so-called single molecular magnets^, the energy level diagram consists of discrete levels 
because the molecules contain only small number of magnetic ions and hence the quantum 
dynamics plays important roles. In particular, in the easy-axis large spin molecules such as 
Mn 12 and Fe 8 , step- wise magnetization processes have been found and they are attributed 
to the adiabatic change, that is the quantum tunneling at the avoided level- crossing points, 
and are called resonant tunneling phenomena^. The Landau-Zener mechanism also causes 
various interesting magnetization loops in field cycling processes^. 

The amount of the change of the magnetization at a step in the magnetization process is 
governed by the Landau-Zener mechanism and depends significantly on the energy gap at 
the crossing. This dependence has played an important role in the study of single-molecule 
magnets. Observations of the gap have been done on isolated magnetic molecules^. 

The quantum dynamics of systems of strongly interacting systems which show quantum 
phase transitions is also of much contemporary interest. As far as static properties are 
concerned, the action in the path-integral representation of a <i-dimensional quantum system 
maps onto the partition function of a (d + l)-dimensional classical model, which is the key 
ingredient of the quantum Monte Carlo simulation. 8 From this mapping, it follows that the 
critical properties of the ground state of the <i-dimensional quantum system are the same as 
those of equilibrium state of the (d + l)-dimensional classical model, quantum fluctuations 
playing the role of the thermal fluctuations at finite temperatures. 

However, from a view point of dynamics, the nature of the quantum and thermal fluc- 
tuations are not necessarily the same. Thus, it is of interest to study dynamical aspects of 
quantum critical phenomena. As a typical model showing quantum critical phenomena, in 
the present work we adopt the one-dimensional transverse Ising modeR 

Recently, interesting properties of molecular chains which are modeled by the transverse 
Ising model with large spins have been reported^. However, in this paper, we concentrate 
ourselves in systems of S = 1/2. The dynamics of the transverse Ising model plays important 
roles in the study of the quantum annealing in which the quantum fluctuations due to the 
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transverse field are used to survey the ground state in a complex system^. The dynamics of 
domain growth under the sweep of the transverse field through the critical point has been 
studied related to the Kibble-Zurek mechanism^^. 

In this paper, we study the hysteresis behavior as a function of the external field in 
the ordered state by performing simulations of pure quantum dynamics, that is by solving 
the time-dependent Schrodinger equation. 14 This gives us numerically exact results of the 
dynamical magnetization process of the transverse Ising model under sweeps of magnetic 
field with constant velocities. 

Previously we have studied time evolution of magnetization of the transverse Ising model 
from a view point of Landau-Zener transition, sweeping the field slowly and finding transi- 
tions at each avoided level crossing point.- However, for fast sweeps the transition at zero 
field H z = disappears and the magnetization does not change even after the field reverses. 
The magnetization remains in the direction opposite to the external field for a while, and 
when the magnetic field reaches an certain value, the magnetization suddenly changes to 
the direction of the field. This sudden change is also found for very slow sweeps at the level 
crossing point. However, the present case has the following two differences: (1) the switching 
field does not necessarily corresponds to a level crossing, and (2) in all cases the changes 
are independent of the size L of the system. This sudden change resembles the change of 
magnetization at the coercive field in the hysteresis loop of ferromagnetic systems, where it 
is called spinodal decomposition. Therefore, we will call the phenomenon that we observe in 
the quantum system a "quantum spinodal decomposition" and the field "quantum spinodal 
point" Hsp. We study the dependence of Hsp on the transverse field H x , and also study the 
sweep- velocity dependence of H SP . 

As in the cases of the single molecular magnets, it should be possible to understand the 
dynamics of the magnetization in terms of the energy levels as a function of field. However, 
because the structure of the energy-level diagram strongly depends on the size of the system, 
it is difficult to explain the size- independent property of the quantum spinodal decomposition 
from the energy-level structure only. In the case of much faster sweeps, we find almost 
perfect size-independent magnetization processes. We also find a peculiar dependence of 
magnetization on the field in the case of weak transverse fields. These processes can be 
understood from the energy-level diagram for local flips of spins, but not from the energy 
diagram of the total system. 
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In this paper, we attempt to understand the microscopic mechanism that gives rise to this 
size independent dynamics. We introduce a perturbation scheme for fast sweeps, regarding 
the fast sweeping field term as the unperturbed system and treating the interaction term 
as the perturbation. From this viewpoint, we investigate fundamental, spatially local time- 
evolutions which yield the size-independent response to the sweep procedure. In particular, 
we propose a perturbation scheme in terms of independent Landau-Zener systems, each of 
which consists of a spin in a transverse and sweeping field. A system consisting of locally 
interacting Landau-Zener systems explains well the magnetization dynamics for fast sweeps. 



II. MODEL 

We study characteristics of dynamics of the order parameter of the one-dimensional 
transverse- Ising model with periodic boundary condition under a sweeping field. 9 The Hamil- 
tonian of the system is given by 

W(t) = - J £ o?o* +1 -H X J2 of - H z (t) Y: v?, (1) 

i i i 

where of and of are the x and z components of the Pauli matrix, respectively. Hereafter, we 
take J as a unit of the energy. The order parameter is the z component of the magnetization 

M* = £<. (2) 

i 

We study dynamics of the order parameter of the model, i.e., the time dependence of the 
magnetization under the time dependent field H z (t) 

(M z ) = (*(t)|M z |*(t)), (3) 
where is a time dependent wavefunction given by the Schrodinger equation 

tt^|*(t)> = W(f)|tt(t)>. (4) 

In the present paper we study the case of linear sweep of the field 

H z {t) = -H + ct, (5) 

where — Hq is an initial magnetic field. In the present paper, we set H /J= 1, and c is the 
speed of the sweep. We use a unit where % — 1. 



In the case H z = 0, the model shows an order-disorder phase transition as a function of 
H x . The transition point is given by H£ = J. In the ordered phase (H x < J), the system 
has a spontaneous magnetization m s : 

m s = lin l lim(G(0)|M*|G(0)>, (6) 

H z — >+0 L— >oo 

where |G(0)) is the ground state of the model with H z = 0. Therefore, the ground state is 
twofold degenerate with symmetry-broken magnetization, while the ground state is unique 
when H x > J. Because of these twofold symmetry-broken ground states, the magnetization 
changes discontinuously at H z = 0. 

In a finite system L < oo, this degeneracy is resolved by the quantum mixing (tunneling 
effect) and a small gap opens at H z = 0. This gap becomes small exponentially with L 
as shown in Appendix A. Therefore, the change of the magnetization becomes sharper as 
L increases. Dynamical realization of this change by field sweeping becomes increasingly 
difficult with L. This phenomenon corresponds to the existence of metastable state. 

The energy-level diagram becomes complicated when L increase. However, as shown 
below, when c is large the system shows a size-independent magnetization dynamics which 
is not easily understood in terms of the energy-level diagram. In this paper, we focus on the 
regime of moderate to large sweep velocities. 



III. ENERGY STRUCTURE 



In Fig. deleft), we present an energy- level diagram for L = 6 and H x = 0.7. We plot all 
energy levels as a function of H z . We find that the energy levels show a linear dependence 
at large fields, where quantum fluctuations due to H x have little effect. The levels are mixed 
in the region —3 < H z /J < 3 where the energy levels come close and are mixed by the 
transverse field. The isolated two lowest energy levels are located under a densely mixed 
area, which represent the ordered states with M = L and — L, and they cross at H z = with 
a small gap AE\, reflecting the tunneling between the symmetry broken states. The gap 
AE\ is so small that one cannot see it in Fig. [Deleft). After the crossing, these states join 
the densely mixed area. In Fig. [T^right), we show the energy levels for L = 16 where we plot 
only energies of a few low-energy states. In this figure, we also find the above mentioned 
characteristic structure of two lowest energy levels. 
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Let us point out a few more characteristic features of the energy-level diagram. A finite 
gap AE 2 exists between the crossing point of the low-lying lines and the densely populated 
region of excited states. The d- dimensional Ising model in a transverse field is closely related 
to the transfer matrix of (d + l)-dimensional Ising model. From this analogy, we associate 
AEi to symmetry breaking phenomena. When symmetry breaking takes place, the two 
largest eigenvalues of the transfer matrix of the model become almost degenerate. The 
energy gap corresponds to the tunneling through the free-energy barrier between the two 
ordered states and vanishes exponentially with the system size. On the other hand, AE 2 
is related to the correlation length of the fluctuation of anti-parallel domains in the order 
state. The correlation length is finite at a given temperature in the ordered state and is 
almost size-independent. At H z = we can calculate eigenenergies analytically, and we can 
explicitly confirm that AE± vanishes exponentially with L and that AE 2 is almost constant 
as a function of the size. The dependencies of the energy gaps at H z = are discussed in 
Appendix A. 

For large H z , the slopes of the low lying isolated lines are ±L because they represent 
the states with M = ±L. Thus, the field at which the lines merge in the area of densely 
populated excited states is given by 

AEo 

h z ^^ = h;(l). (?) 

At this point, the magnetization shows a jump when the speed of the sweep is very slow 4 
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However, as we will see in the following sections, the dynamical magnetization does not 
show any significant change at this field value when the sweeping field is fast. Another type 
of jump will occur that we called quantum spinodal jump or quantum spinodal transition. 

IV. EVOLUTION OF THE MAGNETIZATION FOR FAST SWEEPS OF THE 
FIELD 

A. Quantum spinodal decomposition 

When we sweep the magnetic field from H z — — 1 to H z — 1, the magnetization shows 
a rapid increase to a positive value. In Fig. [2j we depict examples of dynamics of the 
magnetization as a function of time for a sweeping velocity c = 0.001. Because H z (t) = 
—H + ct, H z also represents time. 

The magnetization stays at a negative value until a certain field strength is reached. The 
system can be regarded as being in a metastable state. Then, the magnetization changes 
significantly towards the direction of the field in a single continuous jump, the magnetization 
processes M z {t) depending very weakly on the system size. In the classical ordered state, 
we know a similar behavior. Namely, at the coercive field (at the edge of the hysteresis), 
the magnetization relaxes very fast and the relaxation time does not depend on the size. 
Thus, we may make an analogy to the spinodal decomposition phenomena. We call the 
phenomenon that we observe in the quantum system "quantum spinodal decomposition" 
and we call the field at which the magnetization changes Hsp- It should be noted that the 
spinodal decomposition corresponds to the fact that the size of the critical nuclei becomes 
of the order one. If the size of the critical nuclei is larger than the size of the particle as 
in the case of nanoparticles, the critical field of the sudden magnetization reversal, which is 
also a kind of spinodal decomposition, strongly depends on the size. 

First, let us attempt to understand this dynamics from the view point of the energy 
diagram. As we mentioned in the previous section, the low-lying levels of M = ±L merge 
with the continuum at H z = H*. Thus, we expect that at this point the magnetization 
begins to change because the states with M = ±L begin to cross other states. In fact, in 
earlier work, we found stepwise magnetization processes at avoided level-crossings in very 
slow sweeps, each of which could be analyzed in terms of successive Landau-Zener crossings. 3 
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FIG. 2: (color online) Left: Magnetization M z {t) as a function of H z (t) for c = 0.001, H x = 0.5 
and various system sizes. Solid (red) line: L = 12; Dashed (green) line: L = 14; Dotted (magenta) 
line: L = 16; Right: Same as left except that H x = 0.7. 



From Fig.[2peft), we find that the sharp change of M z (t) starts at H z = 0.2 ~ 0.25, which 
is much larger than H*(L). We estimate H*(12) ~ 0.18, and for larger lattices H*(L) is even 
smaller for larger lattices. Moreover, it should be noted that the magnetization processes 
display almost no size-dependence. In Fig. fright), which shows M z (t) for H x = 0.7, we also 
find that the magnetization processes M z (t) for all sizes L are very similar. Here H$p ~ 0.11 
is again significantly larger than H*(L) (for L — 16 and AE 2 ~ 1.4 in the case H x = 0.7, 
and hence if* (14) ~ 0.09). This observation is in conflict with the picture based on the 
structure of the energy level diagram given earlier. 

In Fig. [31 we present an example of sweep-velocity dependence for a system with L = 20 
(results of other sizes are not shown). The magnetization processes show strong dependence 
on the sweep velocity c, as expected. However, for fixed c, there is little dependence on L 
(results of other sizes are not shown). 

We have found the characteristic change in the cases of relatively large quantum fluc- 
tuations, i.e. H x = 0.5 and 0.7. The size-independence indicates that the change occurs 
locally. When H x is small, the quantum fluctuations are weak and local flips of clusters 
consisting of small number of spins become dominant. In Fig. HJ M z (t) for H x = 0.1 is 
shown, where a peculiar sequence of jumps is found. It is almost independent of the system 
size (except for L = 2). Before the large jump of the magnetization at H z / J — 1, there is a 



9 




-1 -0.8-0.6-0.4-0.2 0.2 0.4 0.6 0.8 1 

H z /J 

FIG. 3: Magnetization M z (t) as a function of H z (t) for L = 20, H x = 0.7 and various sweep 
velocities. Solid (red) line: c = 0.1; Dashed (green) line: c = 0.01; Dotted (magenta) line: 
c = 0.001. 

small but non-zero precursor jump around H z ~ 2/3. After these jumps, the magnetization 
shows a plateau of M z (t)/L w — 1/2 until the smooth crossover to the saturated value takes 
place around H z /J = 2. The value H z /J = 2 corresponds to the spinodal point of the 
corresponding classical model. 

The positions of these jumps can be understood from the viewpoint of local " cluster" flips. 

Let us consider a single spin flip, that is, a flip from the state with all spins | ■•) 

to a state | 1 •••). The diabatic energies of these states are E = —LJ + LH Z 

and Ei = — (L — 4) J — (L — 2)H Z , respectively. Thus, the crossing of these states occurs at 
= 4 J/2 = 2 J. The transition probability due to the transverse field H x at this crossing 
is proportional to H^., because the matrix element for a single flip is proportional to H x . 

If we consider a collective flip of a connected cluster of m spins, the diabatic energy of 
this state is 

E rn = -(L-4)J+(L-2m)H 2 , (8) 
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FIG. 4: Magnetization M z (t) as a function of H z (t) for H x = 0.1, c = 0.001, and various system 
sizes. Solid (black) line: L = 2; Solid (red) line: L = 4; Long dashed (green) line: L = 6; Dashed 
(magenta) line: L = 8; Dotted (red) line: L = 10; Dashed dotted (blue) line: L = 12; 



and thus, the crossing of the states occurs at 

H (m) = 4J / 2m = 2J/m. (9) 

For m — 2, 3, ... we have if]" 1 ) = 1, 2/3, . . . respectively. These values do not depend on L. 
It should be noted that for the system L = 2, the 2-spin cluster (m=2) surrounded by + 
spins can not be realized, and no jump appears at H z / J = 1. 

The matrix element for the m-spin cluster flip is proportional to H™ (see Appendix A). 
Therefore, for small H x , only the flips with small values of m are appreciable. In the case of 
H x /J = 0.1 for c = 0.001, jumps for m < 2 are observed. The change of the magnetization 
of each spin is given by a perturbation series and is independent of L as shown in Appendix 
B. These local flips may correspond to the nucleation in classical dynamics in metastable 
state. 

If c becomes small or H x becomes large, contributions from large values of m become 
relevant. Then, magnetization process consists of many jumps, and amount of the change 
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becomes large. But, as long as the perturbation series converges, we have a size-independent 
magnetization process, as shown in Fig. [2j This sharp and size-independent nature is con- 
sistent with the property of the classical spinodal decomposition. 

In the classical system, the magnetization relaxes to its equilibrium value at the spinodal 
decomposition point. In contrast, for pure quantum dynamics, the magnetization of the state 
does not change for adiabatic motion along a particular energy level. Only if we include an 
effect of contact with the thermal bath, relaxation to the ground state takes placed. 



1. Phase diagram 

In Fig. [5j we give a schematic picture of the order parameter M as a function of the tem- 
perature T and the field H z in the thermal phase transition of a ferromagnetic system. The 
overhanging structure signals the existence of the metastable state. The spinodal point is at 
the edge of the metastable branch. In this figure, the magnetic field is swept from positive 
to negative, and the metastable positive magnetization jump down to the equilibrium value 
at H SP (T). 

In a mean field theory for the magnetic phase transition at a finite temperature, the 
spinodal point is given by 



H SP = -Jz\ll- — + —]n ^_ Ji _ !BjS . \, (10) 



j 

where z is the number of nearest neighbor sites. We show the dependence of H$p as a 
function of T by a dash-dot curve in Fig. [5j 

A similar argument can be made for the classical ground state energy. Let the z- 
component of spin be denoted by a. Then, the energy is expressed by 

E = -Ja 2 - H x Vl - a 2 - Ha. (11) 

We assume that the energy satisfies the condition 

dE , s 

which gives 

-2Ja+-=^=-H = 0. (13) 
V 1 — cr 
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FIG. 5: (color online) Schematic picture of the magnetization M as a function of field H below 
the critical temperature. Open circles denotes the spinodal decomposition points. The (red) dash- 
dotted curve in the H-T plane shows Hsp(T) as given by Eq. (jTUJ) . 

Here, we consider the metastable state and thus we set H = —\H\ for a > 0. At the end 
point of metastability, 

da dH , 

_ = o or — = 0. (14) 



This leads to 



tt 2/3\ 

' 1 1 - 1 If ) ■ < 15 > 



The end point of the metastable state is given by 

/ , tt N 2/3\ 3 /2 

which gives H$p as a function of H x and is shown in Fig. O as the long-dashed curve. 

It is interesting to note that expression Eq. ffTBT) is very similar to the well know expression 
of the Stoner-Wohfarth model^ for the reversal of a classical magnetic moment under the 
application of a magnetic field tilted with respect of the easy anisotropy axis. This is not 
surprising because, with both a longitudinal and transverse field component, this model can 
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be considered as a realization of the classical spinodal transition. One might derive the 
Stoner-Wohfarth model from equation Eq. ffTTl) by replacing the exchange energy parameter 
J by the anisotropy energy constant D. 

It should be noted that the critical H x is a factor of two larger than that of the correct 
value = J for the one- dimensional quantum model. This difference is due to the presence 
of quantum fluctuations. Therefore, in Fig. [6] we plot Eq. (fl6]) with and without renormalized 
values of the fields. The long-dashed curve denotes the case of H x scaled by 1/2, and the 
dashed curve denotes the case where both H x and H z are scaled by 1/2. 

As we saw in Fig. [2J we find a large change of magnetization at a values of H z for each 
value of H x , which we called i?sp. In Fig. 0, we plot values of H z at which (1) M(t) shows 
a small but clear jump, (2) M(t)/L is equal to —1/2, and (3) M(t) saturates as a function 
of H z , for various values of c. The data show a dependence on H x that shows a similar 
dependence to the dotted line. If we use other value of c, the values of H z change. Although 
the values of H z for (1), (2) and (3) for larger values of c are larger than those for c = 0.001, 
the values of H z for c = 0.0001 are close to those for c = 0.001. They seem to saturate around 
the value of the dotted line, and we may identify a sudden appearance of size independent 
change as an indication for a quantum spinodal point. If we sweep much faster, the jumps 
of the magnetization becomes less clear, as we now study in more detail. 

B. Very fast sweeps 

For a fast sweep c = 0.1, the magnetization processes for different sizes almost overlap 
each other, see Fig. [7J The data for L = 14, 16, and 20 are hard to distinguish. This 
almost perfect overlap is rather surprising from the viewpoint of the structure of energy- 
level diagram. The data for L = 6 deviates from the others. This fact indicates that for 
these parameters (H x = 0.7, c = 0.1) the relevant size of the cluster (m in Eq. ([8])) is larger 
than 6 but smaller than 14. 

Let us now study the behavior if we sweep much faster. In Fig. [Speft), we show the 
magnetization as a function of t (or H z (t)) for L = 12 with c = 10, 20, 50, 100 and 200. For 
these parameters, the data for other L are almost indistinguishable from the L — 12 data 
and are therefore not shown. As Fig. [8](left) shows, the magnetization oscillates about a 
stationary value for large values of H z where the energy levels with different magnetization 
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FIG. 6: (color online) Spinodal points H~! as a function of the quantum fluctuation H x for various 
sweep velocities c. The horizontal dotted lines correspond to 2/m for (to = 2, . . . , 10) (see Eq. ([9])). 
Plusses (1), crosses (2), and stars (3): c = 0.0001; Open squares (1), solid squares (2), and solid 
diamonds (3): c = 0.001; Open circles (1), bullets (2), and open diamonds (3): c = 0.01; Open 
triangles (1), solid triangles (2), and inverted solid triangles (3): c = 0.1. The numbers (1), (2), 
and (3), correspond to the field at which M(t) shows a small but clear jump, M(t)/L = —1/2, and 
M(t) saturates as a function of H z , respectively. 



M are far separated in the energy-level diagram as we saw in Fig. Q3 Let us study the c- 
dependence of the saturated value Ms = Ms(c). In Fig. fright), we plot the change of the 
magnetization AM/L = (Mg(c) — (—L))/L as a function of 1/c. As shown in Fig. [fright), 
the data can be fitted well by the expression 

A*A*. + . (17) 

±j ±j c 

where AMq / L and a are constants. These constants, to good approximation, do not depend 
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FIG. 7: (color online) The magnetization M(t) as a function of the sweeping field H z for H x = 0.7, 
c = 0.1, and various system sizes. Solid (black) line: L = 6; Solid (red) line: L = 14; Long dashed 
(green) line: L = 16; Dashed (magenta) line: L = 18; Dotted (red) line: L = 20; Dashed dotted 
(blue) line: L = 12. Except for L = 6, all other curves overlap, indicating that for sufficiently large 
systems, the dependence on L is very weak. 

on the system size. 

In order to explain the observed 1/c dependence, we introduce a perturbation scheme for 
fast sweeps (see Appendix B). We regard the sweeping field (Zeeman) term as the zero-th 
order system and treat the interaction among spins as the perturbation term. The result is 
a series expansion in terms of H /c (see Eq. (IB8I) ). which explains the 1/c dependence. 

In Appendix B, we also introduce a perturbation scheme based on independent Landau- 
Zener systems each of which is given by a spin in a transverse field with a sweeping field. 
In Appendix B, we show that this scheme can explain the behavior of the magnetization 
dynamics in the fast-sweep regime. 
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FIG. 8: Left: The magnetization M z (t) as a function H z for L = 12, H x = 0.7 and various sweep 
velocities. Solid (red) line: c = 10; Long dashed (green) line: c = 20; Dashed (magenta) line: 
c = 50; Dotted (dark blue) line: c = 100; Dashed dotted (magenta) line: c = 200. Right: AM/L 
as a function of 1/c for Hx = 0.7 and various system sizes. Solid (red) line: L = 4; Long dashed 
(green) line: L = 6; Dashed (dark blue) line: L = 10; Dotted (magenta) line: L = 12. 



V. SUMMARY AND DISCUSSION 

We have studied time-evolution of the magnetization in the ordered phase of the trans- 
verse Ising model under sweeping field H z . We found significant jumps of the magnetization 
at a certain value of the magnetic field which we called quantum spinodal point H z p . Al- 
though the energy-level diagram of the system significantly changes with the system size, 
we found size-independent magnetization processes for each pair (H x ,c). 

In principle, it should be possible to understand the quantum dynamics of the magneti- 
zation from the energy-level diagram of the total system. Indeed the picture of successive 
Landau-Zener scattering processes works in slow sweeping case^. However, for fast sweeps, 
the time evolution can be regarded as an assembly of local processes, the interaction between 
the spins being a perturbation. Hence the dynamics of the magnetization does not depend 
on the size. 

When the quantum fluctuations are weak (small H x ), a series of local spin flips governs 
the magnetization dynamics. The jumps of magnetization can be understood on the basis 
of energy- level crossings of certain spin clusters (Eq. ([8])). The energy- level structure cor- 
responding to the local cluster flips is, of course, present in the energy-level diagram of the 
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total system but it is hidden in the complicated structure of the huge number of energy 
levels. 

For large values of H x and fast sweeps, the magnetization process is also size-independent. 
To explain this feature, we have introduced a perturbation scheme in which the small pa- 
rameter is Hq/c. In addition, we introduced a new perturbation scheme based of single-spin 
free Landau-Zener processes, which all together have allowed us to provide an understanding 
of the magnetization dynamics under field sweeps in terms of the energy-level structure. 
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APPENDIX A: SIZE DEPENDENCE OF THE ENERGY GAP AT H z = 

The eigenvalues of the model Eq. (TTJ are given by-2 

£7 = (2^-1), (Al) 

g 

where rj and ry are fermion annihilation and creation operators, respectively, and 



c^ = 2^1 + 2Acosg + A 2 , (A2) 
where A = J/H x . When the number of the fermions is even, q takes the values 

and when the number of the fermions is odd 

2tt 4tt tt(L-2) /k . 

q = 0, ±— , ±— , • ■ ■ , ± \ \ 7T- (A4) 

Because u q > 0, the ground state is given by 77^77 = 0. Thus, in the case of even number of 
fermions, the ground state is given by 



L/2 



E El = -2H X ]T x j 1 + 2A cos ( L + A 2 , ( A5) 



m=l 
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and the first excited state is given by 



Kv;> = E m : l#, v /:i - 2Acus(^^) A 2 . (AG) 



In the case of an odd number of fermions, the lowest energy state is 



m=l 

and the first excited state is given by 



E i = H x Vl + 2A + A 2 + H x Vl - 2A + A 2 

-2H X 1 ^/l + 2Acos(^^^') + A 2 , ( A7) 



L -2 



E02 = Eoi + 2H X V1 - 2A + A 2 + y 1 + 2A cos ( — —j + A 2 . (A8) 
The energy gaps are given by 

AE l = E Q1 -E m , (A9) 

AE 2 = E E2 - E E1 , (A10) 



and 



AE 3 = E Q 2 - E E1 . (All) 



Using these formulae, we can calculate the L-dependence of the gaps. The results are plotted 
in Fig. [9](left). We find that AEi vanishes exponentially with L, that is 

AEi oc exp(— aL), (A12) 

where a depends on A. We also plot — logA£'i/2 to confirm the exponential dependence. 
On the other hand, we find that AE 2 is almost independent of L, and AE 3 is very close 
to AE 2 , reflecting the fact that above the 3rd level the infinite system has a continuous 
spectrum. 

It is also of interest to study the dependence of the energy gap AE\ on H x for several L. 
In Fig. fright) we show the data on double logarithmic scale. In the regime of small H x we 
find a linear dependence on H x , suggesting that 

AEx oc Hf. (A13) 

Indeed, for small H x the slopes of the lines is given by 2S' = L. This dependence on H x and 
L is to be expected when L spins flip simultaneously. 
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FIG. 9: (color online) Left: Size-dependence of the energy gaps for H x = 0.7. Bullets: Difference 
500 x AEi/J between the energy of the first excited state and the ground state energy. This 
difference vanishes exponentially with L. Solid squares: — (In A£i/J)/2; Stars: AEijJ\ Open 
squares: AE3/J. Right: The energy gap AEi/J as a function of H x for several L. Plusses: 
L = 10. Crosses: L = 20; Stars: L = 30; Open squares: L = 40. Note the double logarithmic 
scale. In both figures, lines are guides to the eyes only. 

APPENDIX B: PERTURBATION ANALYSIS FOR LANDAU-ZENER TYPE 
SWEEPING PROCESSES 

When the sweep velocity c is very large, the duration of the sweep is very short. This 
suggests that it may be useful to study the magnetization processes by a perturbational 
method in terms of the small parameter 1/c. 

Let us consider the following model. 

H = H + ctV, (Bl) 

where Tio and V are time independent. We will work in the interaction representation with 
respect to ctV, that is, we take the motion of ctV as reference, not Ho as is usual done. The 
Schrodinger equation is 

ih%-\*) = (?io + ctV)\*). (B2) 
at 

In the interaction representation we have 

I*) = e - ict2v/2h \$), (B3) 
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and the equation of motion is given by 

= ih(-ictV/K)e~ lct2y/2h \§) + ihe- ict2v/2h ^\<$>) = e ~ ict2v / 2n (ctV + ih^J |$>, 

(B4) 

and therefore the Schrodinger equation for |<3>) is given by 

ih%-\<!>) = e tct2v ^ n H e- tct2v / 2h m. (B5) 

Defining 

W{t) = e ict2v / 2h H e- ict2v/2h , (B6) 
we can use the usual perturbation expansion scheme for 

ih^\*)=W(t)\*), (B7) 



and find 



!*(*)> 



1 + Q i w ^ dt + Q 2 i £ + 



|$(0)>. (B8) 



In the sweep (—Ho < ct < H ), to = —H /c and t = H /c. Thus, the integral is of order 
Ho/c. Therefore we can regard the above expansion is a series expansion in terms of power 
of Hq/c. Of course the series can be also regarded as a power of Ho as in the usual sense. 

1. Transverse Ising model under a field sweep 

Now, we consider our problem 

H(t) = - J £ a]a] +1 -H x J2<?--ctJ2 <>]■ (B9) 

3 3 3 

We set 

Wo = -JE^+i-^E^ ( B1 °) 

3 3 

and 

v = -E°s- ( Bn ) 

Then, W(t) is given by 

= - J E^M + i "^E + .-e^ 2 /") . (B12) 
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We may include the diagonal term —JJ2j a j a j+i m V- Then the expansion is regarded as 
series of H x . This expansion corresponds to the series of jumps discussed in Eq. (jSJ). 

We also note that if J = the above process is an ensemble of independent Landau-Zener 
processes. Each of them is independently expressed by 



inJ^I$ LZ (*)> 



-H x ( 



(B13) 



2. Perturbation theory in terms of independent LZ systems 



Next, we consider the case in which the transverse field is included in V. We sweep the 
field from — H to H . The duration of the sweep is 2H Q /c. We assume that 



#o > J> H x , 



(B14) 



such that the motion due to V is that of an ensemble of independent Landau-Zener processes. 
Thus, we consider the ensemble of the LZ systems as the unperturbed system. 

We know the properties of each system. Namely, we know that the scattering becomes 
small when c becomes large. The time evolution of each LZ system is given by 



f 1 ] 


/ 






\ 




(B15) 



in the adiabatic basis, that is in the representation that uses the eigenstates of the system 
with given H z (t). Here, p is the probability for staying the ground state. In the Landau- 
Zener theory, p is given by the well-known expression 

( 7i m\ 



p 



cxp 



he 



(B16) 



In the case of small Hq, p may have a different form. Even in those cases, the expression 
Eq. (1B15I) is still correct and the present formulation works if we employ a correct expression 
for p. 



The unperturbed state is given by 



M*)=n^)= e 



iL<f>{t) 




<g> • • • ® 



^/T=pe i£ ® 

fair) 
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FIG. 10: Energy level diagram for a two-spin Landau-Zener model Eq. (|B18p with H x = 0.7. Left: 
J = 0; Right: J = 1. 



The zero-th order is given a usual Landau-Zener process of which the energy diagram is 
given by Fig. fTOT left). which shows the energy- level diagram for the two independent spins. 

Thus, in this case there are four states, (++), (H — ), ( — h), and ( ). The states consisting 

of (H — ) and ( — h) are degenerate with energy zero. 

The interaction term — JJ2j v'j&j+i is the perturbation. As far as the expansion Eq. (1B8[) 
converges with less than L-th terms, it gives a local effect. To first order in J, only the 
nearest-neighbor spins interact, giving a contribution of the order J. The sweep- velocity 
dependence is taken into account through the zero-th order term. If we take a large H , 
the integral in Eq. flB8j) is no longer small, and we have to regard Eq. flB8j) as a series of 
J. Therefore, we do not have any small parameter, and the Eq. (1B8I) represents the original 
general dynamics. In the case of fast sweeps, the effective range of quantum mixing in which 
the diabatic energy levels (levels for H x = 0) cross each other, is of order Lx J, and therefore 
the duration of interaction is of order LJ/c. Hence, the integration gives a contribution of 
order LJ/c which now becomes the small parameter. In the case of finite H , the small 
parameter is the minimum of (H Q /c, LJ/c). In the present study, H = 1. Then H /c is the 
small parameter and we cannot use the form of p given in Eq. flB16j) . In any case, the series 
converges for the fast sweeps and we expect that the perturbation effect does not depend 
on L. 

The system described by the first-order perturbation theory corresponds to a Hamiltonian 
of two spins exhibiting the Landau-Zener scattering process and which are coupled by an 
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Ising interaction. The Hamiltonian reads 

ftcLZ = -J°l°l ~ (H x af - cta z x ) - (H x a x 2 - cta z 2 ) . (B18) 

The energy-level diagram of this system is shown in Fig. [TUT right). Let us study the effect 
of the interaction on the dynamics in this case. We compare the magnetization processes of 
the model Eq. (IB18I) with J = and J = 1. The results are shown in Fig. [TTT left) . Note 
that the sweep starts from H z = —H = — 1. 

Next, in Fig. [TTlfright). we show the magnetization processes for c = 100 for the model 
Eq. (IB18I) with that of the same model with J replaced by 2 J. If we use a small value of 
H Q , the ground states of the models at H z = — 1 differ significantly. Therefore, to compare 
the results, in this figure, we take H = —60 such that the ground state of both models is 
close to the all-spins-down state. The average of the first and the third curves is close to 
the second curve. This fact indicates that the processes are well described by the first-order 
perturbation theory. Indeed, the deviation from the single Landau-Zener model is 0, J, and 
2J, respectively. 

We also compare the magnetization processes of the models Eq. flB18j) with J replaced 



by 2 J and that of a model with 3 spins in Fig. [121(a). The difference between the models 
of 4 spins and of 12 spins is also shown in Fig. [T27 right). In all these cases, we start at 
H z = — 1 because the magnetizations per spin are very close in all the cases. We find almost 
no difference, indicating that the processes are well described by the first-order perturbation 
theory. 

When the sweep velocity becomes small, we may need higher order perturbation terms. 
If the relevant order of the perturbation is less than the length of the chain, we expect a size- 
independent magnetization process. The size independent magnetization in the quantum 
spinodal decomposition can be understood in this way. 

The local motion of magnetization can be understood from a view point of an effective 
field from the neighboring spins. We may study the magnetization process of a single-spin 
in a dynamical mean-field generated by its neighbors.— Let us describe the situation by the 
following Hamitonian: 

H MF = -(H z (t) + 2J(a z })a z - H x a x . (B19) 

Because the mean field is almost 2 J during the fast sweep, the mean field simply shifts H 
by a constant 2J. Thus, we conclude that for fast sweeps, the dynamics is very similar to 
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FIG. 11: Left: Comparison of the magnetization processes of the model Eq. (jB18[) with J = 
(thin line) and J = 1 (thick lines) for H x = 0.7. Solid (red) line: c = 10; Long dashed (green) line: 
c = 20; Dashed (magenta) line: c = 50; Dotted (dark blue) line: c = 100; Dashed dotted (magenta) 
line: c = 200. For each c, the magnetization of the single LZ process is shifted by an amount such 
that at H z = — 1 it coincides with the magnetization of the model Eq. (|B18P with J replaced by 
2 J. Right: Comparison of the magnetization of a single LZ process, that of the model Eq. (|B18|) . 
and that of the model Eq. ([518]) with J replaced by 2J. H x = 0.7 and c = 100. Solid (black) line: 
c = 10, single LZ process; Solid (red) line: c = 10, Eq. (|B18P : Long dashed (green) line: c = 10, 
Eq. (|B18|) with J replaced by 2 J; Dashed (magenta) line: c = 100, single LZ process; Dotted (red) 
line: c = 100, Eq. (jBl~8~]) : Dashed dotted (blue) line: c = 100, Eq. (lBl"8i) with J replaced by 2 J. 



that of a single spin, meaning that for the dynamics, the effective field on each spin in the 
lattice is essentially that same as the applied field. This conclusion is consistent with our 
earlier comparison of the zero-th and first-order perturbation results. 
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